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GENERATION OF A DESIRED WAVEFRONT WITH A PLURALITY OF PHASE 
CONTRAST FILTERS PHASE 

FIELD OF THE INVENTION 

The present Invention relates to a method and a system for synthesizing a prescribed 
intensity pattern based on generalized phase contrast Imaging. 
BACKGROUND OF THE INVENTION 

It Is well known to form an image by phase contrast Imaging methods In which phase 
modulation of light i s converted into Intensity modulation. As opposed to intensity 
modulation, phase modulation does not involve loss of energy. 
A generalized phase contrast imaging method and system is disclosed in 
international Patent Application Publication No. WO ge/34207 for synthesizing a 
prescribed intensity pattern. The generalized method is not based on the so-celled 
Zemike approximation that the phase shift * is less than 1 radian. An Improved 
method is provided without this assumption and based on an imaging operation with 
a s.mple one-to-one mapping between resolution elements or pixels of a spatial 
phase modulator and resolution elements of the generated intensity pattern. 

HHZZS PhaM ,ma9in9 meth ° d * s ^*"9 a prescribed Intensity 

pattern l(x',y*), comprises the steps of: 

• dividing the intensity pattern l(x-./) into pixels In accordance with the disposition 
of resolution elements <x.y) of a spatial phase mask having 

a plurality of Individual resolution elements (x,y). each resolution element (x,y) 
modulating the phase of electromagnetic radiation incident upon it with a 
predetermined phasor value e^, 

• radiating electromagnetic radiation towards the spatial phase mask. 

• Fourier or Fresnel transforming the modulated electromagnetic radiation, 

• Phase shifting with a spatial phase filter (SPF) in a region of spatial frequencies 
comprising DC in the Fourier or Fresnel plane, the modulated electromagnetic 
radiation by a predetermined phase shift value 9 in relation to the remaining part 
of the electromagnetic radiation, and 

• forming the Intensity pattern by Fourier or Fresnel transforming, respectively the 
Phase shifted Fourier or Fresnel transformed modulated electromagnetic 
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radiation, whereby each resolution element (x,y) of the phase mask is imaged on 
a corresponding resolution element (x'.y) of me image, 

• calculating the phasor values «r*» of (he phase mask and the phase shift value 
9 in accordance with 

for selected phase shift values 6. a being the average of the phasora e^-" of 
the resolution elements of the phase mask. 

• selecting, for each resolution element, one of two phasor values which represent 
a particular grey level, and 

• supplying the selected phasor values e**« to the resolution elements (x,y) of the 
spatial phase mask. 

In one embodiment disclosed in WO 98/34207, the spatial phase mask is positioned 
at the front focal plane of a lens while the spatial phase filter Is positioned In the back 
focal plane of the lens, whereby a first electromagnetic field at the spatial phase 
mask is Fourier transformed by the lens Into a second electromagnetic field at the 
spatial phase Alter. Thus, specific spatial frequencies of the first electromagnetic field 
w,ll be transmitted through the spatial phase filter at specific positions of the phase 
filter. For instance, the energy of the electromagnetic radiation at zero frequency 
(DC) ,s modified by the phase filter and transformed onto the intersecting point of the 
Founer plane and the optical axis of the lens also denoted the zero-order diffraction 
region. 

Typically, the spatial phase filter is adapted to phase shift the DC-part of the 
electromagnetic radiation and to leave the remaining part of the electromagnetic 
radiation unchanged. Alternatively. theDC-part of the electromagnetic radiation 
remains unchanged and the remaining part of the electromagnetic radiation is phase 
shifted. This alternative is preferred when the eneigy level of the DC-part of the 
electromagnetic radiation Is so high that the phase shifting part of the phase fitter will 
be destroyed by it For example in laser cutting, the DC level of the laser beam can 
be so high that a phase shifting dot positioned at the intersecting point of the DC part 
of the laser beam at the phase filter would deteriorate. It is also possible to block the 
electromagnetic radiation (no transmlttance) In the zero-order diffraction region, 
however, the DC energy of the radiation Is then lost 
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SUMMARY OF THE INVENTION 

The above-mentioned method operates on a plane Incoming electromagnetic field. It 
Is an object of the present invention to provide a method and a system operating on a 
more complex electromagnetic field for further flexibility, e.g M to be able to utilize light 
5 sources that are not based on a single coilimated point source, such as laser arrays, 
VCSEL-arrays, LED-arrays etc. 

Further, the spatial phase filter according to the known method processes an incident 
electromagnetic field in the very small diffraction region only, e.g. the zero-order 
diffraction region, leaving the remaining part of the filler un-utilized. Distributing the 
1 0 spatial processing of an electromagnetic field across a larger area of the spatial 

phase filter would make the filter capable of handling higher energy levels. Further, it 
should be possible to match the phase shifting regions of the spatial phase filter to a 
source for emission of electromagnetic radiation that cannot be considered to 
constitute a point source, because it emits radiation from a large area. 

15 According to a first aspect of the present Invention, the above and other objects are 
fulfilled by a phase contrast system for synthesizing an intensity pattern i(x\ y) f 
comprising 

• a source of electromagnetic fields for emission of at least two substantially plane 
electromagnetic fields with different axes of propagation. 

20 • a phase modifying element for phase modulation of the electromagnetic fields by 
phasor values e*^ and positioned so that the at least two electromagnetic 
fields are incident upon it at different respective angles of approach, 

• first Fourier or Fresnel optics for Fourier or Fresnel transforming the phase 
modulated electromagnetic fields positioned in the propagation paths of the at 

25 least two phase modulated fields.' 

• a spatial phase filter with at least two phase shifting regions positioned at 
respective zero-order diffraction regions of the at least two respective phase 
modulated electromagnetic fields for individually phase shitting the at least two 
respective Fourier or Fresnel transformed electromagnetic fields by 

30 predetermined respective phase shift values e n In relation to the remaining part 

of the at least two respective transformed electromagnetic fields, and 

• second Fourier or Fresnel optics for forming the intensity pattern y') by 

4 Fourier or Fresnel transforming the at least two respective phase shifted Fourier 

or Fresnel transformed electromagnetic fields. 
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the phasor values e**** of the phase modifying element and the phase shift 
values 0 n substantially fulfilling that 

for selected phase shift values 0 n , wherein 

5 A is an optional amplitude modulation of the spatial phase filter outside the zero- 
order diffraction regions, 

B« is an optional amplitude modulation of the spatial phase filter In the respective 
n'th zero-order diffraction region, 

a = (a|exp(z$k) te the average of the phasors e*™ of the resolution elements 
10 of the phase modifying element, and 

^ = and 

S(n) is the Intensity of the n'th electromagnetic field, and 
K n = I - J 0 (1 .22*% ) , wherein 
Jo Is the zero-order Basse! function and 
15 7i„ relates the radius R 1n of the n'th zero-order filtering region to the radius R 2 of 

the main-lobe of the Airy function of the input aperture, r\ n « ^ iR* . 

According to a second aspect of the present invention, the above and other objects 
are fulfilled by a method for synthesizing an intensity pattern l(x\ /), comprising the 
steps of 

20 • dividing the Intensity pattern Kx'Y) Into pixels In accordance with the disposition 
of resolution elements (x.y) of a phase modifying element having 

a plurality of individual resolution elements (x,y), each resolution element (x,y) 
modulating the phase of electromagnetic radiation incident upon it with a 
predetermined phasor value 

25 m radiating at least two substantially plane electromagnetic fields with different 
axes of propagation towards the phase modifying element so that the at least 
two electromagnetic fields are incident upon it at different respective angles of 
approach. 
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• Fourier or Fresnel transforming the phase modulated electromagnetic fields. 

• phase shifting In at least two phase shifting regions positioned at respective 
zero-order diffraction regions of the at least two respective phase modulated 
electromagnetic fields for Individually phase shifting the at least two respective 

5 Fourier or Fresnel transformed electromagnetic fields by predetermined 

respective phase shift values e n in relation to the remaining part of the at least 
two respective transformed electromagnetic fields, and 

• forming the intensity pattern K*. y') by Fourier or Fresnel transforming the at 
least two respective phase shifted Fourier or Fresnel transformed 

10 electromagnetic fields, 

• calculating the phasor values e 1 *** of the phase modifying element and the 
phase shift values 0 n substantially in accordance with 

for selected phase shift values 0 m wherein 

15 A is an optional amplitude modulation of the spatial phase filter outside the 

zero-order diffraction regions, 

Bn is an optional amplitude modulation of the spatial phase filter in the 
respective n'th zero-order diffraction region, 

a = |a| exp ) is the average of the phasors e'* x "> of the resolution 
20 elements of the phase modifying element, and 

0 and 

» 

S(n) is the Intensity of the n'th electromagnetic field, and 
K * m 1 - (l^2^7 n ). wherein 
Jo Is the zero-order Bessel function and 
S tf n relates the radius R lfl of the n'th zero-order filtering region to the radius R 2 

of the main-lobe of the Airy function of the input aperture. ti„ = . 

• selecting, for each resolution element, one of two phasor values which represent 
a particular grey level, and 
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• supplying the selected phasor values e^ y> to the respective resolution elements 
(x, y) of the phase modifying element. 

The mathematical expressions will be further explained below. 

The axis of propagation of a plane electromagnetic field is perpendicular to the 
5 electric and magnetic fields. 

Preferably, the at least two substantially plane electromagnetic fields are non- 
interfering Le. the fields are mutualiy incoherent and/or orthogonally polarized and/or 
have different wavelengths, etc. 

It should be noted that, in each resolution element of the phase modifying element, 
10 one of two phasor values which represent a particular grey level of the intensity 
pattern i(x\/) may be selected. 

In an embodiment of the present invention, the spatial phase filter substantially does 
not attenuate the electromagnetic fields incident upon it outside the phase shifting 
regions, i. e. A is equal to one or approximately equal to one. 

15 In an embodiment of the present invention, the spatial phase filter substantially does 
not attenuate the electromagnetic fields incident upon it inside the phase shffting 
regions, I.e. B n Is equal to one or approximately equal to one. 

It is also preferred that the phase, shift value 0„ substantially fulfils the equation 



2|sin^ | 
2 

20 for a loss less filter with A - 1 and 2?, « 1 Jn a preferred embodiment of the present 
invention, the phase shift e n is equal to ir or approximately equal to n. Accordingly 
the previous equation leads to Kjz = y 2 and the phase values &{ x > y) of the phase 
modifying element may be calculated in accordance with 

A 

25 where A is the Illuminated area of the phase modifying element 

The electromagnetic field or radiation may be of any frequency range of the 
electromagnetic spectrum, I.e. the gamma frequency range, the ultraviolet range, the 
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visible range, the Infrared range, the far infrared range, the X-ray range, the 
microwave range, the HF (high frequency) range, etc. The present invention is also 
applicable to particle radiation, such as electron radiation, neutron radiation, etc 

Preferably, the electromagnetic fields are monochromatic or quasi-monochromatic so 
5 that the energy of the electromagnetic fields is concentrated in a narrow frequency 
bandwidth. Since the intensity pattern is reconstructed by, for each of the incoming 
plane fields, interference of two electromagnetic fields generated by different phase 
shifting of different parts of the incoming field. It is required that the frequency range 
of the emitted electromagnetic field is sufficiently narrow to ensure that the two 

1 0 electromagnetic fields are coherent so that their superposition generates the desired 
intensity pattern, if the frequency range is too broad, the two fields will be incoherent 
and the phase information will be lost as superposition of non-coherent fields results 
in a summation of the intensities of the two fields. It is required that the difference 
between individual delays of electromagnetic fields to be superpositioned is less than 

15 the wavelength of the fields. This is a relaxed requirement that allows the 

electromagnetic fields to be relatively broad-banded. For example in the visible 
range, a Xe-lamp or a Hg-lamp can be used as a light source In a system according 
to the present invention with the advantage compared to a laser light source that 
speckle noise is reduced. The requirements of the spatial coherence of the 

20 electromagnetic fields depend upon the space bandwidth product of the 

corresponding system and how close the required system performance is to the 
theoretically obtainable performance of the system. 

Preferably, the electromagnetic radiation is generated by a coherent source of 
electromagnetic radiation, such as a laser, a semi-conductor laser, a strained multi- 

25 quantum well laser, a vertical cavity surface emitting laser (VCSEL), a maser, a 

phase-locked laser diode array, a light.emitting diode, a pulsed laser, such as a sub- 
picosecond laser, etc, or an array of such sources. However, as already mentioned, 
a high-pressure ere lamp, such as an Hg lamp, a Xe lamp, etc. may also be used and 
even an incandescent lamp may be used as a source of electromagnetic radiation. 

30 The phase modifying element changes the phase of an electromagnetic field Incident 
upon it Optionally, it may also change the amplitude of an electromagnetic field 
Incident upon it. The phase modifying element may transmit or reflect the incident 
electromagnetic field. The phase modifying element is divided into a number of 
resolution elements, each of which modulates the incident electromagnetic field by 

35 changing Its phase by a specific predetermined value. The predetermined values are 
assigned to each resolution element In different ways depending upon the technology 
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applied In the component For example In spatial light modulators, each resolution 
element may be addressed either optically or electrically. The electrical addressing 
technique resembles the addressing technique of solid-state memories In that each 
resolution element can be addressed through electronic circuitry to receive a control 
signal corresponding to the phase change to be generated by the addressed 
resolution element The optical addressing technique addresses each resolution 
element by pointing a light beam on it, the intensity of the light beam conespondlng 
to the phase change to be generated by the resolution element Illuminated by the 
light beam. 

Spatial phase modulation may be realized utilizing fixed phase mask, a liquid crystal 
device based on liquid crystal display technology, a MEMS (micro electro-mechanical 
system), a MOEMS (micro opto-electro-mechanical system), such as a dynamic 
miiror device, a digital micro-mirror array, a deformabie mirror device, etc. a 
membrane spatial light modulator, a laser diode array (integrated light source and 
1 5 phase modulator), smart pixel arrays, etc 

Seiko-Epson produces a transmitting liquid crystal SLM (LCSLM) having a high 
resolution matrix of transparent liquid crystal elements wherein the relative 
permittivity of each element can be electrically modulated In order to vary the 
refractive index and thereby the optical path length of the element. 
Meadowlark produces a parallel-aligned liquid cjystal (PAL-SLM) with a high fill 
factor, but this device has a very low resolution In that It contains only 1 37 phase- 
modulating elements. 

Hamamatsu Photonics produces a dynamically controllable PAL-SLM with VGA or 
XGA resolution. 

Texas Instruments produces a Digital Mirror Device (DMD) having an array of mirrors 
each of which can be tilted between two positions. 

In an embodiment of the invention, the phasor values are contained In an object that 
is not controlled by the system, e.g. a specimen In a slide, in stead they are detected 
In the system. 

The spatial phase filter is typically a fixed phase mask, such as an optically flat glass 
plate coated with a dielectric layer In the region wherein the modulated 
electromagnetic field is phase shifted 6 n in relation to the remaining part of the 
electromagnetic field. However, the spatial phase modulators mentioned In the 
previous section may also be used for spatial phase filters. In addition, non-linear 



20 



25 



30 



28/08 2003 TOR 15: 38 FAX +45 33970071 ALB I HNS A/S 



©010/033 



6 



10 



15 



20 



26 



30 



materia,* providing for self-phase modulat/on such as KenMype materials can a(so De 
used for introducing phase shift <?„. 

An Imaging system maps the phase modulating resolution elements of the phase 

Tit? ° n ^ "** SUrfaCe ° f ^ reconstructed '"tensity pattern. It may 

comprise a 4f-lens configuration (two Fourier transforming lenses utilizing 
transmission of light or one Fourier transforming tens utilizing reflection of light) or a 
angle .maglng lens. However, any optical imaging system providing a filtering plane 
for the spatial phase filter may be applied in a phase contrast imaging system. 
In the > method and system according to the present invention, the intensity pattern 
1 I ,S , Senerated by su P e T**ltlon of electromagnetic fields in the image plana of 
the Imaging system. The phase modifying element changes the phase values of an 
electromagnetic field incident upon It and the Imaging system directs the 
electromagnetic field with changed phases reflected from or transmitted through the 
Phase modifying element towards the spatial phase niter. The phase filter phase 
shifts a part of the electromagnetic field and the imaging system is adapted to 
superimpose in the Image plane the phase shifted part of the electromagnetic field 
wfth the part of the electromagnetic field that Is not phase shifted by the spatial phase 
Tiiier* 

feZ^n H TT 6d emb0d,ment * inVenti ° n ' *• Ph8Se »«**« — « 

.s positioned at the front focal plane of a lens while the spatial phase filter is 
positioned in the back focal piane of the lens, whereby a first electromagnetic field at 
the phase modrfying element is Fourier transformed by the lens into a second 
electromagnetic field at the phase fitter. Thus, specific spatial frequencies of the first 
electromagnetic field wfll be transmitted through the spafial phase filter at specific 
positions of the phase filter. For Instance, the enemy of the electromagnetic field at 
zero frequency (DC) Is transmitted through the phase filter at the Intersecting point of 
theFouner plane and the optical axis of the lens also denoted the zero-order 
diffraction region. 

is an important advantage of the invention that utilisation of arrays of sources is 
facilitated in that the positioning and/or shapes of the phase shifting regions of the 

ITvTZT* be *° *• 9e ° mel,y ° f the ** ""ear 

array of VGSELs forms the source, the phase shifting regions of the spatial phase 

filterform a corresponding linear army of phase shifting regions, each of the regions 

being positioned at the zero-order diffraction region of a respective VCSEL in the 
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VCSEL array. Further, the shape of each phase shifting region may match the shape 
of the zero-order diffraction region of the respective VCSEL. 

Likewise, a phase filter may match a source with a specific geometrical shape with a 
continuous phase shifting region covering an area of the phase filter that corresponds 
5 to the zero-order diffraction region of the source. 

Thus, the energy of the electromagnetic fields of the system may be distributed over 
a large area compared to the area of a zero-order diffraction region of a single plane 
electromagnetic field of a known phase contrast imaging system. 
Thus, the phase shifting regions of the spatial phase filter may form a rectangular 
array, a circular array, a linear array, two linear crossing arrays, a continuous region, 
a ring, etc. a 

The at least two substantially plane electromagnetic fields with different axes of 
propagation may be generated in a time multiplexed manner. e.g. by a scanning 
m.rror or prism, deflecting or reflecting a beam of electromagnetic field in different 
1 S directions of propagation. 

The capability of handling high energy levels of electromagnetic fields of the present 
invention may be utilized for provision of a laser cutter comprising a system 
according to the present invention. 

Further, the capability of handling high energy levels in combination with the 
capability of generating crossing light beams may be utilized for provision of an 
optical tweezer or an array of optical tweezers according to the present Invention, 
in an embodiment of the present invention, wherein the apertures of the system is of 
Insignificant Importance to the operation of the system and calculation of the phasor 
values, i.e. K„ is equal to one or approximately equal to one. 
For a more detailed understanding of the Invention, the Zemike approximation is 
reviewed below, followed by a generalization where the above-mentioned 
mathematical expressions axe derived for a single (n=1) on^xls centred phase 
contrast filtering implementation. 

The Zemike phase contrast method allows for the visualization of phase 
perturbations by the use of a Fourier plane phase shifting filter. The Dutch physicist 
Fritz Zernike received the Nobel Prize in 1953 for inventing this method, which led to 
a break-through in medicine and biology by making essentially transparent cell or 
bacteria samples clearly visible under a microscope. Its successful operation 
however, requires that the spatial phase distribution, tf*.,), at the input is limited to 
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a "small-scale" phase approximation where the largest phase is typically taken to be 
significantly less than */3. If the phase distribution at the Input is thus restricted, 
then a Taylor expansion to first order te sufficient for the mathematical treatment so 
that the Input wavefront can be written as 

5 «p(4(*»?)H+ (1) 

The light corresponding to the two terms in this "small-scale" phase approximation 
can be separated spatially by use of a single lens where the phase distribution is 
located in the front focal plane and the corresponding spatial Fourier transformation 
is generated in the back focal plane of the lens. In this first order approximation the 

10 constant term represents the amplitude of on-axfs light focused by the lens in the 
bacK focal plane and the second spatially varying term represents the off-axis light. 
Zemike realized that a small phase shifting quarter wave plate acting on the focused 
light makes it possible to obtain an approximately linear visualization of small phase 
structures by generating interference between the two-phase quadrature terms in Eq. 

15 (1): 

It should be noted that a three-quarter waveplate works equally well to produce 
contrast but the plus sign in Eq. (2) is negated leading to so-called negative phase 
contrast. A substantial improvement in the visibility of the Zemike phase contrast 
visualization in Eq. (2) requires strong damping of the focused light in addition to the 
phase shift required to generate the contrast. 

In the general case, where we are not limited to a small-scale input phase 
perturbation we cannot assume that a series expansion to first order as in the 
Zemike approximation is a sufficient representation of a given phase perturbation. 
Higher order terms In the expansion need to be taken into account, so the expansion 
takes the form: 

However, here the spatially varying terms can not be considered as separate from 
the supposedly focused light represented by the first term In this Taylor series 
30 expansion, as Is Implied by the Zemike approach, and all of these spatially varying 
terms contribute to the Intensity of the on-axis focused light. For a significant 
modulation in the input phase, this contribution of the spatially varying terms can 
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result In a significant modulation of the focal spot amplitude In the back focal plane of 
the lens. These terms can in fact result in either constructive or destructive 
Interference with the on-axis light, although the net result will be an attenuation of the 
focused light amplitude, which only has a maximum value for a perfect unperturbed 
plane wave at the input. 

For phase objects breaking thB Zernlke approximation we must, therefore, find an 
alternative mathematical approach to that of the Taylor expansion given in Eq. (3). 
We have chosen a Fourier analysis as a more suitable technique for completely 
separating the on-axis and higher spatial frequency components. This gives the 
following form for exp(i^ x ,y)) t where (x,y) eCl: 



«P(*(**))- [fl** /^(^(^^J^+^gherfiequcncytcrnis" (4) 



in this Fourier decomposition the first term is a complex valued constant linked to the 
on-axis focused light from a phase object defined within the spatial region, Q , and 
the second term describes light scattered by spatially varying structures In the phase 
object Comparing Eq. <3) and Eq. (4) It is apparent that the first term of Eq. (3) is a 
poor approximation to the first term of Eq. (4) when operating beyond the Zemike 
small-scale phase regime. 

An Important issue to consider when analysing the effect of spatial filtering on the 
light diffracted by phase perturbations is the definition of what spatially constitutes 
focused and scattered light. In the previous description of Zemike phase contrast it 
was assumed that the focused light is spatially confined to a somewhat unphyslcal 
delta function. As we know, any aperture truncation such as occurs within a practical 
optical system will lead to a corresponding spatial broadening of the focused light It 
is therefore essential that we define the terms "focused light" and "scattered light- 
explicitly for such a system. In this context it is necessary to look more carefully at 
the sequence of apertures confining the light wave propagation through a typical 
optical set-up. 

A commonly applied architecture that provides an efficient platform for spatial filtering 
is illustrated in Fig. 1 and is based on the so-called 4-f configuration. An output 
interferogram of an unknown phase object or phase disturbance Is obtained by 
applying a truncated on-axis filtering operation In the spatial frequency domain 
between two Fourier transforming lenses (L1 and t2). The first lens performs a 
spatial Fourier transform so that directly propagated light Is focused Into the on-axis 
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fllterlng region whereas spatially varying object information generates light scattered 
to locations outside this central region. We can describe a general Fourier filter in 
which different phase shifts and amplitude damping factors are applied to the 
"focused" and "scattered" light. In Fig 1. we show a circulariy symmetric Fourier filter 
5 described by the amplitude transmission factors A and B for the "scattered" and 
-focused" light respectively and by the relative phase shift ft These filter parameters 
can be chosen to replicate any one of a large number of commonly used filter types 
(Le. phase contrast, dark central ground, point diffraction and field-absorption 
filtering). By applying a given Fourier filter and a second Fourier lens, we obtain an 
10 interference pattern in the observation plane. The focused on-axis light acts as the 
synthetic reference wave (SRW) in the CP! system, this interferes with the scattered 
light to generate the output interference pattern. In the following section we discuss 
the importance of the SRW and show how it influences, among other things, the 
choice of the Fourier niter parameters. 

15 Having described the generic optical system that makes up the CPI. we turn to a 
detailed analytical treatment of the Important elements in this system. Assuming a 
circular input aperture with radius. Ar . truncating the phase disturbance modulated 
onto a colllmated. unit amplitude, monochromatic field of wavelength, X . we can 
describe the Incoming light amplitude a(x,y) by, 

20 «(^)=circ(//Ar)exp(^(^,^)) (5) 

at the entrance plane of the optical system shown in Fig. 1 using the definition that 
the circ-function Is unity within the region, r « v^7"i Ar . and zero elsewhere. 
Similarly, we assume a circular on-axis centred spatial filter of the form: 

H {f»fyh « l P(^)-l)oirc(/ r /V r )] (6) 

25 where B e[0; 1] is the chosen filter transmlttance of the focused light. 0 «[<*2x] is 
the applied phase shift to the focused light and A € [0; l] is a filter parameter 
describing field transmittance for off-axis scattered light as indicated In Fig. 1. The 
spatial frequency coordinates are related to spatial coordinates in the filter plane 
such that: (/„/,) - Wr(x„y,) and f r ^fifTf*. 

.30 Performing an optical Fourier transform of the Input field from Eq. (5) followed by a 
multiplication with the filter parameters in Eq. (6) and a second optical Fourier 
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transform (corresponding to an inverse Fourier transform with Inverted coordinates) 
we obtain an expression for the intensify /(V.y) describing the interferogram at the 
observation plane of the 4-f set-up: 



10 It should be noted that to achieve a tractable analytic expression In Eq. (7) it has 
been assumed that the spatial frequency content of the phase object is sufficiently 
described by the term, a . within the on-axls centred filtering region characterized by 
the spatial frequency range 6f r , 

The generally complex valued and object dependent term, a , corresponding to the 
1S amplitude of the focused light plays a significant role in the expression for the 
interference pattern described by Eq. (7). Referring to the discussion In the 
introduction, we are now able to confirm that the frequent assumption, that the 
amplitude of the focused light is approximately equal to the first term of the Taylor 
expansion in Eq. (1), can generally result in misleading interpretations of the 
20 interferograms generated at the CPI output. 

Of similar importance in the analysis of Eq. (7) is the term g(r') describing the spatial 
profile of the SRW, diffracted from the aperture formed by the on-axis centred filtering 
region. It Is the interference between this SRW term, carrying the Information about 
the filtering parameters, and the Imaged phase object that generates the output 
25 Interferogram. Thus. H is important to obtain an accurate description for the SRW and 
thereby an accurate derivation for Eq. (7). The zero-order Hankel transform followed 
by a series expansion in the spatial dimension, f, will be used to describe the SRW. 
This is a relatively simple approach, which to the knowledge of the author has not 
previously been applied to this problem. 

30 For a circular Input aperture with radius, Ar , we can describe the radius of the 

corresponding central phase shifting region of the Fourier filter (characterized by the 
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parameters B and in terms of a radial spatial frequency range A/, . We can thus 
obtain the following expression for the SRW by use of the zero-order Hankel 
transform: 

*(r')»2nArf' l y 1 (2rtAj/^ 0 (2nr'/ r )d/ r (9) 

5 In order to simplify the analysis, w© Introduce a term n , whteh explicitly relates the 
radius of the central filtering region, *, , to the radius of the main-lobe of the Airy 
function, R, , resulting from the Fourier transform of the circular Input aperture alone. 
We can thus express rj in terms of Ar and A/ r such that; 

r) = R,/R 2 = (0.6 1)' 1 ArA/, fl 0 ) 

10 where the factor of 0.61 arises from the radial distance to the first zero crossing of 
the Airy function corresponding to half of the Airy malntobe factor, of 1 .22. If we make 
this substitution in Eq.(9) and then perform a series expansion in r' , we obtain the 
following expression for the SRW: 

g(r') = 1 -/„ (l.22«n) [(0.611X71)* A (l .22*n)](r '/Ar) a + 

{[(0.61*7,)' / 4 ][2/, (i.22mi)-o.6imt/, (i.22*7,)]}(r'/Ar) 4 .... < 11 > 

In this expansion, the SRW is expressed in radial coordinates normalised to the 
radius of the imaged input aperture. This can easily be scaled to allow for a 
magnification within the imaging system, though for the remainder of our analysis a 
direct imaging operation is assumed. From Eq.(H) it is apparent that the SRW will 
change as a function of the radius of the central filtering region. Additionally, it is 
clear that the SRW profile is not necessarily flat over the system output aperture. 
This is an Important, yet often neglected, factor in determining the performance of a 
CPI. 

Depending on the accuracy needed for the description of the Interferograms one can 
choose to include a number of spatial higher order terms from the expansion in Eq. 
(11). The influence of the higher order terms has the largest Impact along the 
boundaries of the imaged aperture. For r, -values smaller than 0.627 and when 
operating within the central region of the image plane, spatial higher order terms are 
of much less significance and we can safely approximate the synthetic reference 
wave with the first and space invariant term in Eq. (1 1 ): 



15 



20 



25 
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g (r • e central region)* l-J 0 (l .22itT|) (12) 
so that we can simplify Eq. (7) to give: 

/(*',>' , )s^|« P (^$(x , ,>'•))+A:|5|(fi J 4• , ex P (,e)-^)| ^ (13) 

where K= 1-J 9 (l22trrj). The influence of the finite on-axis filtering radius on the 
5 focused light is thus effectively included as an extra -filtering parameter" so that the 
four-parameter filter set (A,B,0 t K( n )) together with the complex object dependent 
term, a. effectively defines the type of filtering scheme we are applying. 
Having determined a suitable operating range for the CR in terms of the production 
of a good SRW. we must now examine the role that the remaining filter parameters 
10 play in the optimisation of a CPI. 

From Eq. (13) we see that the filter parameters (A,B, B) can be combined to form a 
single complex valued term. C, the combined filter term, such that: 

C = |C|exp (f Vc ) - BA-* exp (»9) - 1 ( 14) 
therefore, Eq. (13) can be simplified to give: 

1S ^^-^KW'^O-^O+irHKf d5) 



where 



.+"±(*rJ*Wrin(y e ) (16) 

20 Since it is a complex variable, the combined filter term C. which effectively describes 
the complex filter space, can be considered to consist of a vector of phase Vc and 
length |Cj as shown in Eq.(i4). Thus in order to obtain an oveiview of the operating 
space covered by all the possible combinations of three Independent filter 
parameters (A,B,e) we can now instead choose to consider a given filter In terms of 

■25 the two combined parameters y, c and |C|. However, referring to Eq.(1 5), it can be 
seen that the filter parameter. A also appears independently of the combined filter 
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term, C. Fortunately, this issue can be resolved by considering that the term BA 
from Eq. (14) must be constrained In the following way: 



BA ml <\ => A = l j? = |c + i| 

BA" »i => w = l. 5 = 1 (17) 



These constraints arise from the adoption of a maximum Irradiance criterion 
minimising unnecessary absorption of light in the Fourier filter, which reduces both 
Irradiance and the signal to noise ratio in the CPI output. 

In the previous sections we derived expressions and graphical analysis methods 
10 relating the spatial average value of a given phase disturbance to obtain peak 
Irradiance and optimal visibility in combination with high accuracy in systems with 
unknown wavefront phase disturbances. We saw that if a CPI is applied to wavefront 
sensing or the visualisation of unknown phase objects the Generalised Phase 
Contrast (GPC) method specifies the filter phase and aperture size parameters for 
15 achieving optimal performance in extracting and displaying the phase information 
carried by the incoming wavefront. On the other hand, in cases where we have 
control over the incoming wavefront or phase modulation the GPC method provides 
extra means of optimisation by encoding the phase distribution itself in addition to 
modifying the filter parameters. The two main scenarios: A) synthesizing the spatial 
phase for intensity display or B) measuring the spatial phase with high accuracy, 
strongly influences which of the parameters in the analysis that should be kept fixed 
and which could be changed or adapted. The first approach is particularly useful 
when the filter parameters have a restricted dynamic range or are fixed. The rigorous 
derivation of the equations for choosing these parameters wlfl be derived in this 
25 section. 

When synthesizing an Input phase distribution for optimal visibility of an output 
intensity pattern the situation is more relaxed than the situation involving accurate 
interferometric measurements of unknown phase disturbances. The parameter n can 
therefore In most cases be chosen to completely encompass the zero-order light with 
the result that the term, K . tends to unity as the Bessel function tends to zero In Eq. 
(12). For this particular case, the SRW becomes a flat top profile where we can 
achieve nearly 100% light efficiency. For smaller and irregular phase patterns fine- 
tuning of v in the region 0.4-0.6 provides for an efficient operating regime while 
maintaining minimal tosses. 
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In order to optimize a synthesized light distribution we wish to generate a dark 
background for the illumination at the observation plane as our principal design 
criterion, so we can write: 

'(*i>^;$o)=o (20) 

5 where indicates observation plane coordinates of the background of the 

generated illumination and ? 0 is the relative phase shift generating a zero-intensity 
level at the observation plane. 

Applying this dark background condition in Eq. (13) we can obtain the following 
expression for a no-loss phase-only filter with filter transmission 
10 parameters. A=*J* = 1: 

*p|(l -exp(ie)) =exp(4) (21 ) 

A key point arising from Eq. (21) is that we now have a simple way of expressing a 
new design criterion relating the spatial average value of any input phase pattern to 
the zero-order phase shift of a matched Fourier phase filter. Since Kis by definition 
1 5 positive and by taking the modulus of Eq. (21 ) we obtain: 

^p|=|2sin(e/2)p {22) 

Eq. (22) is a key result for the fully transmissive wavefront engineered OPC mapping 
thet makes it possible to deduce the range of valid phase parameters fulfilling our 
design criteria from Eq. (20). The largest possible value that the term, tf|a|, takes on 

20 is unity, this leads to the following solution interval for Eq. (22) within a full phase- 
cycle: 

e = [ir/3;5rt/3] (23) 
From Eq, (22) we also observe that can take on a value Jimlted to the interval: 

(24) 

. 25 Eq. (22) and the solution Intervals described by Eqs. (23)-(24) specify the design 

parameters for achieving optimal performance in extracting and displaying the phase 
information carried by the incoming wavefront. Moreover, Eq. (22) hints towards extra 
means of optimisation by encoding the phase modulation depth Itself In addition to 
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the no-toss phase-only filter. This last approach is particularly useful when the filter 
phase has a restricted dynamic range or is fixed. 

Now. assuming that we have a fixed and fully transmlsslve phase-only filter, the best 
choice for the filter parameter is a value that allows for the largest dynamic range of 
phasor values at the input. Accordingly, the smallest possible nee/ value. Koc - 1 / 2 , 
Is desirable implying that 9 * n , leading to the output Intensity distribution: 

/(^y)-2fl-coa(*(* i l y))] (25) 

Inserting Ka « 1 / 2 In Eq. (8) we obtain the following two requirements for the Input 
encoded phase function 4>{x % y) : 

KA' X fcos($(x % y))dxdy = X 

(26) 

A 

We observe that it is only the first requirement in Eq. (26) that Is directly related to the 
output intensity in Eq. (25) via the cosine term. Since there are always two choices 

15 for a given phasor value that result in the same cosine value (excluding 0 and *). we 
notice that the second requirement can subsequently be fulfilled independently of the 
first requirement simply by complex conjugating an appropriate number of phasor 
values. This fact is a key feature of the GPC-method since It makes it possible to 
solely concentrate on the first requirement in the process of synthesizing a desired 

20 an virtually no-loss grey level intensity pattern. 

The first requirement In Eq. (26) can be fulfilled by several means, including: dynamic 
phase range adjustment fill factor encpdlng, phase-histogram adjustment spatial 
scaling of phasor pattern, raster encoding etc In a histogram adjustment technique 
one will typically start out with a desired relative Intensity distribution l{xf 9 y) 4M 

25 where the maximum achievable intensity level is unknown but relative intensity levels 
are known and the lowest intensity level is fixed by the background criterion of Eq. 
(20). The procedure is now to adjust the histogram for l(x',y') <hsvtd while 
maintaining identical relative intensity level ratios until the first requirement in Eq. (26) 
Is fulfilled. Subsequently, the second requirement In Eq. (26) Is fulfilled by complex 

30 conjugating an appropriate part of the phaaors. The simplest procedure fs to complex 
conjugate every second identical phasor value independently of the spatial location. 
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However, this *phasor flipping" procedure can also be turned Into an advantageous 
tool (an extra degree of freedom) for manipulating the spatial frequency content In 
order to optimize the separation of low and high spatial frequency terms at the 
Fourier filter plane by taking the spatial phasor location into account E.g. 
5 neighbouring phasor values can be chosen to have a maximum difference between 
them, thereby Introducing high spatial frequency modulation easing the filtering In the 
spatial Fourier domain. 

In most cases, however, equalized output intensity levels are sufficient In the 
succeeding analysis, we therefore focus on the encoding of the input phase levels to 
1 0 achieve binary output intensity levels, A derivation based on ternary phase levels 
allows for the widest range of binary intensity pattern encoding and automatically 
provides for the simplified but important binary phase level encoding as a special 
case. 

For the ternary phase encoding, we consider the illuminated portion of the Input 
15 aperture area, A, as divided into sub-areas A* A f and A 2 with respective phase 
values # D1 jf, and fa . We are aiming for the derivation of general expressions 
relating the addressing parameters for the phase modulation to the range of possible 
phase parameters of the Fourier filter obeying the design criterion we have already 
set out. We can express the total truncated area and its average phase modulation. 
20 as the sum of the phase-weighted sub-areas: 

4 e*p('$o)+ A cxp (fy )+ A 2 oxp(/$ 2 ) m >4p| (27) 

This can be further simplified by expressing the sub-areas as fractions of the total 
area, A, such that F t e A x I A and F 2 = A 2 /A: 

*i)«P «p(*$,)+ F t exp(^) ^ p| (28> 

25 As previously mentioned we are interested in binary Intensity patterns with levels 
corresponding to the input phase values. In this case the dark background region Is 
defined by (4,,? 0 ) and the bright output level of Intensity, /, Is determined by 

(i4„?,) and at the Input plane. For the binary output Intensity condition it 

follows that: 

30 '(*.)-'(*.) (29) 
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This equality corresponds to a symmetric condition that can be easily verified 
applying the phasor chart analysis technique demonstrated in section 4. Due to this 
symmetry we can simplify the analysis by applying the following substitution: 

A$ = = <ji 0 (30) 

so that Eq.(28) can be rewritten as; 

F % (exp(/A<fr)-l) + R (cxp(-iA*)-l) = K~ } (l -exp(/e))" 1 -1 (31) 

It is now a straightforward task to solve Eq. (31) for the real part and the imaginary 
part respectively, to obtain the following set of equations: 



10 



F t +F 2 =(2K~l)(2K(\-cos(A$))y % 
F t -F 2 =si«(0)(2^sm(A*)(l-cos(e))y 



(32) 



15 



This can also be expressed in terms of the fractional areas, such that 



J* = (AK)' 1 [(2* - - cosCA*))" 1 + gin (8)(sin (A*)(l - cos(6))) M ] 

f 2 »(4i^)" l [(2Ji:•.l)(l-cos(A ^ ^))" , -^(eXsinCA^Xi-cosCe)))" 1 ] 



(33) 



Since we have focused on solutions where identical intensity levels are obtained in 
both the F r region and the F r region we can define the resulting illumination 
compression factor, C, in the following way: 



C=tf +F 2 y l =(i^(2AT l )' , (l-cos(A*)) 



(34) 



20 The minimum compression factor corresponds to uniform illumination at the output 
such that F } + F 2 = 1 , whereas the maximum compression factor Is found to be 
C-*oo for K = y 2 . 

An interesting spedal case can be deduced from Eq. (31) by setting F 2 = 0 . where 
we And that: 



25 



F = F X =(tf(l-exp(ffl))-l)^ (35) 
implying that for the binary phase modulation case we must have: 



28/06 2003 TOR 15:44 FAX +48 39070071 A LB I HNS A/S 



Q 023/033 



-22- 

A<j> = 6 (36) 
In order for the fill factor, F, to be real-valued. 

This result turns out to be the special case that corresponds to the set of solutions 
where a binary phase pattern serves as the Input 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates schematically a known phase contrast Imaging system, 

Fig. 2 illustrates schematically an embodiment of the present invention, 

Fig. 3 also schematically illustrates the embodiment of Fig. 2, and 

Fig. 4 Illustrates various layouts of the spatial phase filter according to the present 
10 invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 shows a known phase contrast imaging system embodied in a 4f common path 
interferometer 1 . A laser emits a light beam, which is expanded by a beam expander 
into a plane light wave of uniform intensity and directs It towards a phase modifying 

1 5 element. The light beam is transmitted through the phase modifying element and a 
Fourier transforming lens (L1). The phase modifying element is positioned in the front 
focal plane of the lens <L1> and a spatial phase filter 6 Is positioned in the back focal 
plane of the tens L1 that is also the front focal plane of a lens L2. The Fourier 
transforming lenses L1, L2 need not have identical focal lengths. Different focal 

20 lengths lead to a magnification ratio different from one. The spatial phase filter 6 
phase shifts by 8. and optionally attenuates (by a factor B), the zero order diffraction 
part 8 of the light modulated by the phase modifying element. Optionally, the 
remaining diffraction part of the light modulated by the phase modifying element may 
be attenuated by a factor A. The reconstructed Intensity pattern l(x\ y') is generated 

25 in the back focal plane 9 of the lens L2. 

Figs. 2 and 3 illustrate schematically an embodiment 10 of the present Invention in a 
4f common path interferometer configuration. An array 12 with a plurality of individual 

light sources PS 1t PS 2 ps„ emits a set of light beams Wave,, Wavez Wave* 

which are colllmated by lens L3 Into a plurality of plane light waves directed towards 
a phase modifying element 14 to be incident upon the element at different respective 
angles of approach. The phase modifying element 14 phase modulates the 
electromagnetic fields by phasor values a!** (x. y) are coordinates of the resolution 
elements or pixels of the phase modifying element 14. The light waves are 



30 
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transmitted through the phase modifying element 14 and a Fourier transforming lens 
Li. The phase modifying element 14 is positioned in the front focal plane of the lens 
L1, and a spatial phase filter 6 is positioned in the back focal plane of the lens L1 that 
is also the front focal plane of a lens L2. The Fourier transforming lenses L1 , L2 need 
5 not have identical focal lengths. Different focal lengths lead to a magnification ratio 
different from one. In each phase shifting region PCF 1t PCF*. ... , PGF n , the spatial 
phase filter 6 phase shifts by respective values 9 U 0* ... . 0* and optionally 
attenuates (by respective factors B 1f B 2 . ... *B«) the respective n'th zero-order 
dftfraction part of the light modulated by the phase modifying element 14. Optionally. 
1 0 the remaining diffraction part of the light modulated by the phase modifying element 
may be attenuated by a factor A. The reconstructed intensity pattern l(x\ y*) is 
generated in the back focal plane 9 of the lens L2. 

The Illustrated system 10 may be controlled by a computer (not shown), comprising 
interface means for addressing each of the resolution elements of the phase 

1 5 modifying element 14 and transmitting a phasor value e 1 ** to the addressed 

resolution element (x. y). Further, the computer may comprise light control means for 
controlling the power of the light-emitting array 12. The computer may also comprise 
input means, such as a keyboard, a diskette drive, an optical disc drive, a network 
interface, a modem, etc, for receiving an image pattern to be synthesized by the 

20 system 10. From the received image pattern, the computer may be adapted to 
calculate phasor values to be transmitted to the resolution elements (x, y) of the 
phase modifying element 14. e.g. based on a histogram technique as described 
herein. Optionally, the phase shifts of the spatial phase filter 6 may be adjustable and 
controllable by optional phase control means of the computer, which may be further 

25 adapted to adjust the phase shift. 

Fig. 4 Illustrates various layouts of the spatial phase filter according to the present 
Invention for illustration of the flexibility provided by the invention. The illustrated 
layouts match corresponding layouts of source arrays. 

PROPOSED APPLICATIONS: 

30 » Laser machine tool e.g. for shaping, forming, assembling, etc, such as cutting, 
drilling, milling, planing, marking, branding, trimming, hardening, scribing, 
labelling, welding, soldering, two- and three-dimensional surfaces especially by 
use of C0 2 and Nd:YAG laser based systems. The main advantage is that 
energy Is not absorbed In the system according to the present invention (thereby 

35 preventing damage of the optical hardware) and this non-absorbed energy Is 
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Instead utilized to Increase the Intensity level of the desired light distribution In 
the Image plane. High power can be delivered to selected regions on a work 
piece simultaneously. 

• Optical tweezer array for manipulation of micro-objects, such as micro- 

5 components, biological cells, etc r using electromagnetic gradient forces 

proportional to the optical intensity pointing in the direction of the intensity 
gradient. 

• Efficient and dynamic spot-array generators based on phase contrast imaging. 
In order to provide bias or holding beams for arrays of optoelectronic elements, 

10 such as bistable elements, photonic switches and smart pixels. 

• Generation of structured light (toss less) for machine vision applications. E.g. 
periodic and skew periodic mesh grid illumination that can be updated In parallel. 

• Photolithographic applications (laser 3D direct writing in parallel without the need 
for sequential scanning). E.g. high power laser direct writing of waveguides in 

15 Ge-doped silica. 

• Spatial light intensity modulation In general by use of pure phase modulation 
(radiation focusators). 

• Laser beam shaping (dynamic). 

• Highly efficient parallel Image projection without the need for a laser-scanning 
20 device. 

• Dynamic Infrared Scene Projection (DIRSP). 

• Exposure device for grating and mask production. 

• LIDAR applications. 

• Laser printing in parallel. 
25 ■ Laser show applications. 

• Atmosphere research. 
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CLAIMS 

1. A phase contrast system for synthesizing an intensity pattern l(x\ y 1 ). comprising 

a source of electromagnetic fields for emission of at least two substantially plane 
electromagnetic fields with different axes of propagation, 

5 a phase modifying element for phase modulation of the electromagnetic fields by 
phasor values ei* M and positioned so that the at least two electromagnetic fields are 
incident upon it at different respective angles of approach, 

first Fourier or Fresnel 

optics for Fourier or Fresnel transforming the phase modulated electromagnetic fields 
10 positioned in the propagation paths of the at least two phase modulated fields, 

a spatial phase filter with at least two phase shifting regions positioned at respective 
zero-order diffraction regions of the at least two respective phase modulated 
electromagnetic fields for individually phase shifting the at least two respective 
Fourier or Fresnel transformed electromagnetic fields by predetermined respective 
1 5 phase shift values 0 n in relation to the remaining part of the at least two respective 
transformed electromagnetic fields, and 

second Fourier or Fresnel optics for forming the intensity pattern l(x\ y') by Fourier or 
Fresnel transforming the at least two respective phase shifted Fourier or Fresnel 
transformed electromagnetic fields, 

20 the phasor values e**^ of the phase modifying element and the phase shift 

values 0„ substantially fulfilling that 

= ^ e*p(«U-l)f 

for selected phase shift values 0 n , wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
25 order diffraction regions, 

B„ Is an optional amplitude modulation of the spatial phase filter In the respective 
n'th zero-order diffraction region, 

a = {ajexp(/^ ) is the average of the phasors e** ( * v) of the resolution elements 
of the phase modifying element, and 



30 



^ = - <f>- , and 
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S(n) is the Intensity of the n'th electromagnetic field, end 
K * - 1 - J* (1 .22*j?„ ) , wherein 
Jo Is the zero-order Bessel function and 

7i„ relates the radius R,„ of the n'th zero-order filtering region to the radius R c of 
the maln-tobe of the Airy function of the Input aperture, t)„ = R, n 

2. A phase contrast system according to claim 1 . wherein the phase shining regions 
of the spatial phase filterform a rectangular array. 

3. A phase contrast system according to claim 1. wherein the phase shifting regions 
of the spatial phase filter form a circular array. 

4. A phase contrast system according to claim 1. wherein the phase shifting regions 
of the spatial phase filter form a linear array. 

5. A phase contrast system according todaim 1, wherein the phase shifting regions 
of the spatial phase filter form two linear crossing arrays. 

6. A phase contrast system according to claim 1. wherein the phase shifting regions 
15 form a continuous region. 

7 A phase contrast system according to claim 6. wherein the phase shifting regions 
of the spatial phase filter form a ring. 

e. A phase contrast system according to claim 1 . wherein the phase shifting regions 
of the spatial phase filter form an arbitraiy array. 



25 



8. A phase contrast system according to any of the previous claims, wherein 



A=1. 



10. A phase contrast system according to any of the previous claims, wherein 



Bn= 1. 



1 1 . A phase contrast system according to any of the previous claims, wherein 



0n a n. 



12. A phase contrast system according to any of the previous claims, wherein 



Ko = 1. 



13. A phase contrast system according to any of the previous claims, wherein the 
source comprises a plurality of light sources. 



© 027/033 



26/06 2003 TOR 15:44 FAX +45 33970071 ALB I HNS A/S 



@ 028/033 



-27- 

14. A phase contrast system according to any of the previous claims, wherein the 
source comprises a laser array, such as a VCSEL array. 

15. A phase contrast system according to any of the previous claims, wherein the 
source comprises a light scanner for time multiplexed emission of the at least two 

5 substantially plane electromagnetic fields with different axes of propagation. 

16. An optical tweezer system according to any of the previous claims. 

17. A laser machining tool according to any of the previous claims. 

18. A phase contrast system according to any of the previous claims, wherein the 
Phaser values are contained in an object that is not controlled by the system. 

10 19. A method for synthesizing an intensity pattern l(x*. /). comprising the steps of 

dividing the intensity pattern Kx'.y 1 ) into pixels in accordance with the disposition of 
resolution elements (x,y) of a phase modifying element having 

a plurality of individual resolution elements (x.y) f each resolution element (x,y) 
modulating the phase of electromagnetic radiation incident upon It with a 
1 6 predetermined phasor value 

radiating at least two substantially plane electromagnetic fields with different axes of 
propagation towards the phase modifying element so that the at least two 
electromagnetic fields are incident upon it at different respective angles of approach, 

Fourier or Fresnel transforming the phase modulated electromagnetic fields, 

20 phase shifting in at (east two phase shifting regions positioned at respective zero- 
order diffraction regions of the at least two respective phase modulated 
electromagnetic fields for individually phase shifting the at teast two respective 
Fourier or Fresnel transformed electromagnetic fields by predetermined respective 
phase shift values 6* in relation to the remaining part of the at least two respective 

25 transformed electromagnetic fields, and 

forming the intensity pattern l(x\ /) by Fourier or Fresnel transforming the at feast 
two respective phase shifted Fourier or Fresnel transformed electromagnetic fields, 

calculating the phasor values e 1 ^ of the phase modifying element and the phase 
shift values 0 ft substantially in accordance with 

30 I(x\y>) s 2 S(»)A> \**p(iHx\y'))+K a «*(«.) -l)f 

for selected phase shift values 0*, wherein 
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A is an optional amplitude modulation of the spatial phase filter outside the zero- 
order diffraction regions, 

B n is an optional amplitude modulation of the spatial phase filter in the respective 
n'th zero-order diffraction region, 

5 a = |5|exp(f^) is the average of the phasors e^> of the resolution elements 

of the phase modifying element, and 
<? = and 

S(n) is the Intensity of the nth electromagnetic field, and 
K n = 1 - J 0 (1 .22**7,, ) . wherein 
10 J 0 is the zero-order Bessel function and 

ij m relates the radius R 1n of the n'th zero-order filtering region to the radius R 2 of 
the main-lobe of the Airy function of the input aperture, r\„ - R 1n I R 2 , 

selecting, for each resolution element, one of two phasor values which represent a 
particular grey level! and 

15 supplying the selected phasor values e'* c *"> to the respective resolution elements (x, 
y) of the phase modifying element. 
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